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Abstract 

Using a sample of 4.7 fb~^ integrated luminosity accumulated with the 
CLEO II detector at the Cornell Electron Storage Ring (CESR), we inves- 
tigate the mass spectrum and resonant structure in t~ — > ir~7r+7r~^'T- de- 
cays. We measure the relative fractions of -ftri(1270) and iri(1400) resonances 
in these decays, as well as the Ki masses and widths. Our fitted Ki reso- 
nances are somewhat broader than previous hadroproduction measurements, 
and in agreement with recent LEP results from tau decay. The larger cen- 
tral value of our measured width supports models which attribute the small 
T~ K~Ti^'K^i'r branching fraction to larger Ki widths than are presently 
tabulated. We also determine the Ka — Kb mixing angle Ok- 
PACS numbers: 13.10+q, 13.35. Dx, 14.40.Aq 
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I. INTRODUCTION 



Decays of the r lepton into three pseudoscalars have been actively studied over the last several 
years. Lately, a number of relatively precise measurements of the branching fractions for t~ — > 
K^TT^TT^u.^ have become available from the ALEPH, CLEO and OPAL collaborations [H-Q. 
However, the resonance substructure of these decays has not yet been measured with high precision. 

The decay r~ — > i^~7r^7r~z^T-, with its simple and well-understood initial state provides in- 
formation on low-Q^ QCD. The effects of SU{3)f symmetry breaking can be observed, the decay 
constants of the Ki resonances can be measured, and the hadronic resonance substructure can be 
studied from an analysis of the final state invariant mass spectra Other interesting topics 

include resonance parameters (such as the widths of the Ki states), tests of isospin relations, and 
measurements of the Wess-Zumino anomaly. Current models of this decay are based on Chiral 
Perturbation Theory (ChPT) calculations. The question of the Ki widths is of special interest 
because the theoretical models based on ChPT provide predictions for the r~ — > K~h^h~i'r 
branching fractions that are significantly larger than current experimental values This dis- 

crepancy can be resolved if the Ki resonances in r decays are much wider than presently measured 
values. In the non-strange sector, it has long been realized that the ai width is considerably larger 
as measured in r~ a^i/j- compared to hadronic production of the oi. The primary goal of this 
analysis is to measure the relative amplitudes of the i^i(1270) and i^'i(1400) resonances that are 
believed to dominate — > K^-K^-K^Ur decays [^,^ and to determine the parameters of the Ki 
resonances. 



II. THEORETICAL ASPECTS OF r ^ iiT tt+vt z^^ DECAYS 

In the Standard Model, the general form for the semileptonic r-decay matrix element can be 
written as 

M = -^n(p,)7'^(l - l5HPr)J^, (1) 

where = {KTnr\V^ — A^\0) is the hadronic current and pu and Pr are the four- momenta of the r 
neutrino and the r lepton, respectively. 

General considerations based on Lorentz invariance and conservation of energy and momentum 
lead to the conclusion that only four independent form-factors are needed to describe the hadronic 
current in r~ K~tt^'k^Vt-. One parameterization I^JH describes this process as 

J'^ = [Fi{si,S2,Q^){p2 - Ps). + F2(S1, S2, Q^){pi - PsUT'^" (2) 
+ Faisi,S2, Q^)e^'''P''piuP2pP3a + S2, Q2)QM, 

where Fi are form-factors, is the Ktttt 4- vector, si is expressed in terms of the final state hadrons' 
momenta pi (i=l for the K~ , i=2 for tt"'" and i=3 for vr") as si = {p2 +^3)^, S2 = {pi +^3)^ and 
T^'^ = {g^'^ — Q^Q'^ /Q"^)- Here, there are two axial vector form-factors Fi and -F2, an anomalous 
vector form-factor Fa, and a scalar form- factor Fg. 



** Charge conjugate states are implied throughout the paper. 
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To derive specific expressions for the form- factors, some assumptions have to be made. It is 
beheved that this decay is dominated by the lowest-mass resonances. There are two axial 
vector resonances which can produce the K~'k~^tt~ final state. These are the weak eigenstates ^Pi 



and ^Pi (us), called Ka and K^- The Kt, couples to the W analagous to a "second class" current, 
violating SU (3) j symmetry. These two weak eigenstates mix with mixing angle 6k to form the 
observable mass eigenstates, Ki(1270) and ii'i(1400) p. The /Ci(1270) subsequently decays into 
K*7T, Kp or Kp', while the i^i(1400) decays almost entirely to K*7r. 
Within the context of ChPT, the form-factors can be written as 

^,^^.BW,„,_,„,SEe^ (3) 

^2 = 7— — ^ BWk*, (4j 

3 (1 + ??) 

where BW denotes a Breit-Wigner mass distribution. The parameter rj is estimated to be 0.33 
in 1^. The coefficients preceding the Breit-Wigner expressions are fixed by ChPT. In the first 
form-factor, the coefficient is taken to be -0.145 based on application of the Conserved Vector 
Current (CVC) to e+e" ^ tt+tt~ data 

In the chiral limit, the scalar form-factor Fg is zero. The vector form-factor Fa is expected to be 
numerically small compared to Fi and F2\ it is only non-zero due to the Wess-Zumino anomaly. The 
vector contribution is approximately 5.5% as calculated using the decay amplitudes found in p. 
For this analysis, we will assume that the vector contribution is zero, and include our uncertainty 
in this term as a systematic error. 

ChPT has found widespread application in tau decays. A model similar to this is used for 
our analysis of the invariant mass distributions in r~ — > K~ 'k'^ 'k~ Vr , as will be described in the 
following sections. 



III. DATA SAMPLE AND EVENT SELECTION 

Our data sample contains approximately 4.3 million r-pairs produced in e'^e~ collisions, corre- 
sponding to an integrated luminosity of 4.7 fb~^. The data were collected with the CLEO II detector 
|TT|| at the Cornell Electron Storage Ring, operating at a center-of-mass energy approximately 10.58 
GeV. 

The CLEO II detector is a general-purpose solenoidal magnet spectrometer and calorimeter. 
The detector was designed for efficient triggering and reconstruction of two-photon, tau-pair, and 
hadronic events. Measurements of charged particle momenta are made with three nested coaxial 
drift chambers consisting of 6, 10, and 51 layers, respectively. These chambers fill the volume from 
r=3 cm to r=l m, with r the radial coordinate relative to the beam {z) axis. This system is very 
efficient (e >98%) for detecting tracks that have transverse momenta {pt) relative to the beam axis 
greater than 200 MeV/c, and that are contained within the good fiducial volume of the drift chamber 
(I cos 6\ <0.94, with 6 defined as the polar angle relative to the beam axis).[^ This system achieves a 



^In this analysis we use charged tracks with momentum above 300 MeV/c. 
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momentum resolution of (Sp/p)'^ = (O.OOlSp)^ + (0.005)^ (p is the momentum, measured in GeV/c). 
Pulse-height measurements in the main drift chamber provide specific ionization (dE/dx) resolution 
of 5.5% for Bhabha events, giving good K/n separation for tracks with momenta up to 700 MeV/c 
and nearly 2a separation in the relativistic rise region above 2 GeV/c. Outside the central tracking 
chambers are plastic scintillation counters, which are used as a fast element in the trigger system 
and also provide particle identification information from time-of-flight measurements. 

Beyond the time-of-flight system is the electromagnetic calorimeter, consisting of 7800 thallium- 
doped Csl crystals. The central "barrel" region of the calorimeter covers about 75% of the solid 
angle and has an energy resolution which is empirically found to follow: 



E is the shower energy in GeV. This parameterization includes effects such as noise, and translates 
to an energy resolution of about 4% at 100 MeV and 1.2% at 5 GeV. Two end-cap regions of the 
crystal calorimeter extend solid angle coverage to about 95% of 47r, although energy resolution 
is not as good as that of the barrel region. The tracking system, time of flight counters, and 
calorimeter are all contained within a superconducting coil operated at 1.5 Tesla. Flux return and 
tracking chambers used for muon detection are located immediately outside the coil and in the two 
end-cap regions. 

We select e"'"e" t^t^ events having a 1-prong vs. 3-prong topology in which one r lepton 
decays into one charged particle (plus possible neutrals), and the other r lepton decays into 3 
charged hadrons (plus possible neutrals). An event is separated into two hemispheres based on 
the measured event thrust axis.^ Loose cuts on ionization measured in the drift chamber, energy 
deposited in the calorimeter and the maximum penetration depth into the muon detector system are 
applied to charged tracks in the signal (3-prong) hemisphere to reject leptons. Backgrounds from 
non-signal r decays and hadronic events with are suppressed by requirements on the impact 
parameters of charged tracks. To reduce the background from two-photon collisions (e^e~ 
e^e~77 with 77 ^hadrons or 77 — > cuts on visible energy {Ems) and total event transverse 

momentum {Pt) are applied: 2.5 GeV < E^s < 10 GeV, and Pt >0.3 GeV/c. We also require the 
invariant mass of the tracks and showers in the 3-prong hemisphere, calculated under the vr^vr^vr^ 
hypothesis, to be less than 1.7 GeV/c. Events are accepted for which the tag hemisphere (1-prong 
side) is consistent with one of the following four decays: r"*" — > e'^i'ei^T, — > fJ-'^t^^J^r, — ^ 7r"'"I7^, 
or r"*" P'^T't- 

Candidate events are distinguished from background r decays with vr'^'s and continuum hadronic 
background (e^e~ qq) by the characteristics of showers in the electromagnetic calorimeter. A 
photon candidate is defined as a shower in the barrel region of the electromagnetic calorimeter with 
energy above 100 MeV having an energy deposition pattern consistent with true photons. It must 
be separated from the closest charged track by at least 30 cm. r~ K~7r~^TT~i'r candidates are 
defined as those events having zero photon candidates in the 3-prong hemisphere. 

The event selection described above provides a sample of events that contains Tr^vr^vr^^'T-, 
T~ K^tt'^tt^Uj- and t~ K~ K^tt~Ut-. In this analysis we neglect possible contributions from 



^The thrust axis of an event is chosen so that the sum of longitudinal (relative to this axis) 
momenta of all charged tracks has a maximum value. 




(5) 
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the decays t~ t:~ K~^tt~u^ and r~ — > K~'k^ K~Uj- because they are unphysical in the Standard 
Model and have not been experimentally observed. We also neglect the r~ K~ K'^ K~Vt- final 
state. This rate is expected to be ~ 1% relative to that for r~ K~'k~^tt~Vt- due to the limited 
phase space and the low probability of forming an (ss) pair from the vacuum. 

IV. RECONSTRUCTION OF INVARIANT MASS SPECTRA 

Due to the very small fraction of kaons in h^h^h^Vj- events]^ and the limited particle iden- 

tification capabilities of the CLEG II detector, it is difficult to identify individual r~ K^tt^tt^Ut 
decays. In this analysis, a statistical approach is used in which the number of t~ K~h^TT~UT- 
events in any given sample is determined using the dE/dx information of the two same-sign tracks 
in the signal hemisphere. The dE/dx analysis is described in detail in [^. 

For each h^h^h^ candidate the invariant mass of the three hadrons is calculated under two 
hypotheses for the first and the third tracks, corresponding to the K~tt^tt~ and tt^tt^K^ mass 
assignments. Each of these two sub-samples is divided into bins of invariant mass. The bins are 
100 MeV/c^ wide, spanning the region 0.8 — 1.7 GeV/c^. After binning in mass, the sub-samples 
that correspond to the same mass bin are combined. 

The dE/dx analysis provides the number of kaons in each mass bin, which is equal to the number 
of r" K~h'^'n~Ur events in that mass interval. The invariant mass spectrum of the K^'k~^'k~ 
system is thereby reconstructed. This distribution contains a contribution from t~ K'^tt^Ut 
decays which must be subtracted, as will be discussed in Sec.^. 

In a similar way the invariant mass spectra of K~'k~^ and vr^vr^ are reconstructed in ten bins 
over the range 0.5 GeV/c^ < Mk-^+ < 1.5 GeV/c^ and 0.2 GeV/c^ < M^+^- < 1.2 GeV/c^, 
respectively. The reconstructed mass spectra are shown in Fig. ||. 

V. BACKGROUND AND EFFICIENCY 

There are two main types of background: continuum hadronic events (e'^'e^ ^ qq ^ hadrons) 
and non-signal r decays. Hadronic background is estimated from a continuum hadronic Monte 
Carlo sample (using the JETSET v7.3 |12] event generator and GEANT [^] detector simulation 



code). This background is subtracted as described in Q. The level of hadronic background is 
approximately 3%. 

r-related background comes primarily from t~ K~ K'^tt~ Uj- decays. These events comprise 
approximately 27% of the events in our reconstructed invariant mass distributions. Smaller back- 
ground contributions arise from t~ — > K~ 'k~ tt^ with incomplete tt^ reconstruction and also 
tau decays to the K~ -k^ tt' ['k^)ut final state through an intermediate K^. These two backgrounds 
comprise 5% and 3% of the events in the K'tt^tt"' invariant mass spectra, respectively. The in- 
variant mass distributions for backgrounds are found using Monte Carlo simulations to obtain the 
shape; the normalization is set by the measured branching fractions P,ll3 • Background predictions 



§Here and later h designates either a kaon or pion. 
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FIG. 1. Reconstructed mass spectra (data) and backgrounds predicted from Monte Carlo sim- 
ulations. 
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are shown in Fig. 0. The invariant mass distributions for ah backgrounds are subtracted from the 
corresponding invariant mass spectra reconstructed from data. 

The efficiency of event reconstruction depends shghtly on the invariant mass. Therefore, it is 
necessary to introduce a mass-dependent efficiency correction. This correction is calculated from r 
Monte Carlo using the KORALB event generator |Q. The maximum variation in efficiency across 
the mass interval of interest is of order 10%. 



VI. FITTING METHOD 

The hadronic structure of the K^tt^tt^ system is investigated by simultaneously fitting three 
invariant mass distributions: M^-^+^-, Mi(-^+ and M^+^-. The fitting function is based upon a 
model similar to the one described in Sec. |, and now outlined in greater detail. 



A. Parameterization of form-factors 

In this analysis, we write the following expression for the axial vector form- factors Fi and F2: 
Fi = 4=(^ • BWu70 + B ■ BWi^oo f^'^^f^"' (6) 

F2 = {-\){C ■ BWi27o + D ■ BWuoo)BWk^ (7) 

that contain the four real parameters A — D. Of these, 3 are independent; the fourth is fixed by the 
normalization requirement that the squared sum of the t~ — > i('^(1270)i/,- and r~ (1400)i'r 
amplitudes must saturate the total t~ K'tt^tt^Vt rate. The coefficients A and B correspond 
to production of the Kp final state through either i('i(1270) {A) or /Ci(1400) (B), modulo a factor 
which includes the appropriate phase space weighting for various final states (denoted as , or 
^^Rb^^)- In our analysis, we fix B to be zero, consistent with current measurements [jl^. Similarly, C 
and D designate production of the K*it final state through the i^i(1270) and -fCi(1400) resonances. 
The decay amplitude parameters in Eqs. therefore correspond to the possible decay chains 

as 

r ^ i^i(1270)i/ ^ Kpv K-ir+TT-Ur : (8) 
r ^ (1270)1/ ^ K*7n^ ^ i^-vr+vr-i/^ : "C" = A ■ • (9) 



r ^ Ki{UOO)iy K*7riy if^vr+vr^i/^ : "L>" = y/l-A"^- (10) 

In Eqs. (P)-(P) we have imposed constraints that follow from the tabulated branching fractions of 
the Ki resonances |l|]: B{Ki{l270) K*tt) = (16 ± 5)% and B{Ki{l270) Kp) = (42 ± 



Thus, in our parameterization of the matrix element one unknown parameter A defines all 
four amplitudes. In addition, the masses and widths of the Ki resonances rj^^(i270)) ^7^^(1400), 
M^j(i270)5 ^i<:i(i400) 3-^6 considered unknown and left as free parameters in the fit. The Breit- 
Wigner distributions for the Ki resonances are defined following the approach of M as 
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BW{s,mK,,rK,) = ^^ ^ ^ 7 . (11) 

The Breit-Wigner distributions for the K* and p resonances contain mass-dependent widths: 

2 

(12) 

where the mass-dependence is defined by Eq. (22) in |^]: 

3 



f P{s) \ ^ 
s \p{ml) 



r(^) = To^ ( ) • (13) 



Here, rriQ and Fq are the nominal mass and width of a particle, p is the momentum of the particle, 
and the variable s is the three- or two-body invariant mass, as appropriate. 

The constants Rx (where X = A,...,D) in Eqs. (P), (P), and (|To|) depend on the masses and 
widths of the 2- and 3-body resonances in this decay and are calculated by numerical integration of 
the appropriate matrix element. Since we are interested in the ratios of quantities (e.g., branching 
fractions), the overall normalization of the Rx parameters is arbitrary. 

In this analysis, we have determined the numerical coefficients for the Breit-Wigner terms 
and -| in Eqs. (|)-(0)] using isospin relations rather than taking the chiral limit as in H. The 
Wess-Zumino anomaly term is set to zero in our model; this term is numerically small enough that 
it can be neglected at our level of accuracy. Appropriate systematic errors are assigned to reflect 
the possible magnitude of this contribution. Note, however, that if the Wess-Zumino anomaly is 
much larger than expected, the K*' may contribute events to the region of M^-^+^- invariant 
mass close to 1.4 GeV/c^, affecting our measurement oi A — D. Note also that we explicitly assume 
all r~ K~Tr^7r~i'r decays proceed through either iri(1270) or Ki{1400). 

In principle, there may be a phase shift between the terms in the form- factors Fi and F2. Such 
phase differences may appear between various decay chains producing the final state K^tt^tt^ . This 
may cause additional constructive or destructive interference and, for example, enhance or suppress 
the p peak in the distribution of vrvr invariant mass. In the most general approach, one would 
introduce three independent phase angles 61, 62 and ^3, corresponding to the possible interfering 
decay chains. However, due to limited statistics, we have neglected such possible interference effects, 
and take into account only the inherent phase of the Breit-Wigner distributions (as described above). 



B. Calculation of the observables 

The interesting observables that we would like to measure are the relative branching fractions 
to the different Ki resonances and the amounts of K* and p in this decay. The decay rate for 
any individual decay chain is proportional to X^Rx- For example, the decay rate for the chain in 
Eq. ^is proportional to A^Ra- In calculating the ratios we choose to normalize to the sum of the 
separate contributions (not including interference effects). With this convention, the fractions of 
different contributions add to 100%. 

With the above definitions we write 

. B{r ^ Ki{mO)iy) ^ A^Ra + C^Rc , 

~ 13{t ^ Ki{1270)iy) + B{t ^ Ki{im)u) A^Ra + C^Rc + D^Rd 
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TABLE I. Results for the full fit. In the fit, the parameter B is set to zero and the parameters C 



and D are constrained by the branching fractions [14| into the K*7r final state from the -ftri(1270) 
and Eri(1400), respectively. The parameters Ra, Rb, Rc and Rd are normalizations returned 
from the fit, as described in the text. Errors for the fit parameters M, T and / are statistical and 
systematic (respectively), as described in the text. 

A = 0.94 ± 0.03 ± 0.04 Ra = 6917 

B = 0.00 (fixed) Rb = 12403 

C = 0.20 (constrained) Rc = 61636 

D = 0.27 (constrained) Rd = 58027 

ri27o = ± 0-08 GeV 

Emoo = O.SOtoil ± 0.14 GeV /1270 = 0.66 ± 0.19 ± 0.13 

M1270 = 1.254 ± 0.033 ± 0.034 GeV/c^ fp = 0.48 ± 0.14 ± 0.10 

Mi4oo = 1.463 ± 0.064 ± 0.068 GeV/c^ 



and 



BiT^KTTTTiy) A^Ra + C^Rc + D'^Rd' 



(15) 



C. Fitting function 

We use a Monte Carlo based fitting procedure, in which a large number (200,000) of simulated 
events are used to simultaneously fit the data distributions for Mk-k-k, ^Kn and M^-tt. From a 
binned fit; we determine the input values of A, Ej^^(i27o)5 rxi(i400)) -^Xi(i270) and Mx^(i4oo) 
which give the best simultaneous match to these mass spectra. As outlined above, our event 
generator is identical to KORALB except that our form-factors [Eqs.(P)-(|^] are used and the 
Wess-Zumino form-factor is set to zero. To take into account finite resolution effects we introduce 
Gaussian smearing of the calculated mass equal to the smearing found from the full GEANT-based 



|13| simulation of the detector. This smearing is typically 5-10 MeV/c . 



VII. FIT RESULTS 

The result of our fit is shown in Figure |2|; the best values for the fit parameters are tabulated 
in Table |. In the same table the values of Rx obtained from numerical integration and the 
derived values for /1270 and fp are also given. The first error in the Table is statistical and the 
second is systematic (discussed in Sec. |Vn| ). The statistical errors on A, the masses and widths 



are calculated using HESSE in MINUIT |]18| and take into account correlations between the fit 
parameters. The asymmetric statistical errors, where appropriate, are also evaluated from the fit. 

The fit results showing contours of constant iii the M1270 vs. M1400 and ri27o vs. ri4oo 
planes are shown in Figures ^ and (Note that the constraint -/Vfi4oo > M1270 is introduced in 
obtaining Figures we therefore fit only over the region above the dashed line.) In these plots. 
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FIG. 2. Fit to the distributions of Mxnn, Mxn and M^^r with individual components over- 
laid, after subtracting the non-signal backgrounds in Fig. 1, and after correcting for the relative 
efficiency as a function of mass. In the upper left plot, the solid line represents the sum of the 
r i<'i(1270)^'r + T ^ Ki{1400)iyr mass spectra and the dashed and dotted lines represent the 
individual i<'i(1270) and i^i(1400) contributions. Similarly, the upper right and bottom plots show 
the K*Tr contributions (dotted) and the Kp contributions (dashed) to the Kmr final state. The 
limited statistics of the fitting function is evident in the non-smoothness of the solid curve in these 
three figures. 12 



3030200-028 




1.0 1.1 1.2 1.3 1.4 

'«K(1270) (GeV/c2) 

FIG. 3. Scan of the mass parameter plane showing fit results for the i^i(1270) mass (horizontal) 
vs. the i^i(1400) mass (vertical). Also shown are 1-4 a standard deviation error contours (statistical 
errors only) around our best fit value (cross) The small open circle shows the present PDG values 
for the iiri(1270) and ifi(1400) masses, with the associated errors (dotted ellipse). The systematic 
errors for our measurement, although not shown, are nevertheless significant (see text). 
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3030200-031 




FIG. 4. Scan of the width parameter plane showing fit results for the i^i(1270) width (horizon- 
tal) vs. the i<'i(1400) width (vertical), as well as contours of constant (1-3 a, statistical errors 
only). The small open circle shows the present PDG values for the Ki(1270) and Ki(1400) widths, 
with the associated errors (dotted ellipse). The dashed line results from the fitting condition that 
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the curves represent la, 2a, 3a and 4(T standard deviation error contours around the best fit point 
(indicated by a cross). 

We perform a secondary fit, in which the masses and widths of the Ki resonances are fixed to 
world average values We obtain /1270 = 0.40 it 0.07, and fp = 0.29 it 0.05 from this second fit 
(statistical errors only), to be contrasted with the significantly larger values extracted from our full 
fit, in which the Ki masses and widths are allowed to float as free parameters. It is not surprising 
that the / values are different in this second fit compared to the original fit, given the high degree 
of correlation between the Ki widths and the relative branching fractions. The for this second 
fit is considerably poorer (30.5/22 degrees of freedom) than the primary fit (12.6/18 degrees of 
freedom) . 



VIII. SYSTEMATIC UNCERTAINTIES 



Systematic errors are summarized in Table 0. The dominant errors are due to the uncertainty 
in the Ki branching fractions and the uncertainty in the background. The errors associated with 
the uncertainty in the branching fractions of the Ki resonances to K*tt and Kp are estimated by 
changing the values for these branching fractions in the form-factors Fi and F2 by one standard 
deviation of the PDG values [Q. The resulting spread in the fit values is taken as the corresponding 
systematic error. 

TABLE II. Summary of systematic error calculations for the full fit. "Ki parameters" refers 
to the uncertainty in the branching fractions, "Model dependence" includes the uncertainty in the 
kaon momentum spectrum in r~ — > K~h~^h~i'r {Pk spectrum). Due to the correlation between 
momentum and three-prong mass, coupled with the fact that we have good particle identification 
only over a limited momentum range, uncertainty in Pk results in a corresponding uncertainty in 
the mass spectrum. Other errors are as indicated. 



Source 




ri270 


ri4oo 




-^1400 


A 


/l270 








GeV 


GeV 


GeV/c2 


GeV/c2 








Ki parameters 




0.005 


0.049 


0.005 


0.004 


0.011 


0.016 


0.012 


Model dependence 


Pk spectrum 


0.057 


0.038 


0.004 


0.020 


0.012 


0.075 


0.055 




Mkktt shape 


0.039 


0.107 


0.024 


0.055 


0.010 


0.083 


0.060 




p' contribution 


0.006 


0.001 


0.005 


0.003 


0.012 


0.005 


0.004 




Vector current 


0.015 


0.026 


0.009 


0.005 


0.005 


0.025 


0.025 


Background level 




0.017 


0.054 


0.020 


0.033 


0.028 


0.054 


0.039 


Function and MC statistics 


0.019 


0.026 


0.005 


0.006 


0.006 


0.026 


0.019 


Total systematic error 


0.076 


0.140 


0.034 


0.068 


0.037 


0.130 


0.097 



The uncertainty in the background is large because of the uncertainty in the branching fraction 
of its largest component, r~ K~ K^'k~Vt-. During the background subtraction the level of all r- 
related backgrounds are varied by amounts corresponding to the errors on the branching fractions 
of these decays [Q. 14]. The hadronic background is similarly varied by 100% to determine the 
systematic error due to our uncertainty in the qq background contribution. 

Another large error comes from the choice of models in our Monte Carlo simulation. This 
includes the uncertainty in the shape of the kaon momentum spectrum {Pk) used to obtain the total 
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number of events with kaons and the uncertainty in the shape of the invariant mass distribution 
for the T~ K~ K^TT^Vr background. These errors are estimated by using several different models 
to extract the invariant mass spectra. For t~ — > K~ tt'^ 'k~ Ur we consider t~ — > Ki{\27{))~Vr 
and r~ Ki{\4Q{))~Vr and the model described in [P; for r~ — > K~ K^tt'i^t, we consider 
r" ^ a5"i/^ ^ K*°K'Ur, p(1690)-i/^ K*^K-Ur and the model described in g. 

There are several fitting function uncertainties. The first is the contribution from the p' which 
may be different in this decay from that observed in e+e" vr+vr" data |1C| due to the phase 
space suppression of p' in our case. Second, the model implemented in our fitting function contains 
no contribution from the vector current. The corresponding fitting function errors from these two 
sources are estimated by varying the level of the Wess-Zumino term and the p' amplitudes from 
zero to the predictions of and [20|. Another possible source of systematic errors is a phase shift 
among the interfering decay chains. In this analysis, the parameters A — D in Eqns. (^)-(|^ are 
real. We have done a study of interference effects with additional phase shifts and found that the 
possible imaginary part of A — -D is consistent with zero at our level of sensitivity. Because the 
fitting function is based upon Monte Carlo, the Monte Carlo statistical error is also included here. 

The bias associated with the procedure of fitting the invariant mass distribution is studied using 
60 samples of signal Monte Carlo with a full detector simulation. The results of this study show no 
systematic shift of the fitted parameter values relative to the input values. Additionally, the errors 
we obtain from analyzing this Monte Carlo sample are fully consistent with statistical expectations. 

In this analysis only the shape of the background-subtracted invariant mass distribution is 
of interest; possible systematic effects that affect the overall normalization of the reconstructed 
spectra are ignored. Among such effects are trigger and tracking efficiencies, and the photon veto. 
Non-r backgrounds (2-photon events, beam-gas interactions, QED background, e.g.) have been 
determined to be negligible for the mass spectrum analysis. 



IX. DISCUSSION OF RESONANCE STRUCTURE 

A. Masses and Widths 

As mentioned previously, theoretical predictions for i3(r — > Kh'^h~Vr) based on ChPT ^ are 
substantially larger than data. However, if the Ki resonances are substantially broader than the 
PDG values, this discrepancy is resolved. In fact, our data suggest larger Ki widths than previous 
world averages (this is evident from Fig. |2|). As indicated in Table ||, we extract the masses 
and widths of the Ki resonances from this fit: Fisro = 0.26l[J:[]? GeV, Fmqo = 0.30+[j:?j' GeV, 
Mi27o = 1.254 ± 0.033 GeV/c^ and Mmqo = 1.463 ± 0.064 GeV/c^. 



In Table III, our result for the widths is compared to the data from ALEPH and DELPHI |15| 
in their analyses of r decays. One observes that all experimental data from r KinTv-r for the 
Ki widths are above the current world averages although the errors remain large. The masses of 
-R'i(1270) and ifi(1400) measured in our analysis are in acceptable agreement with the current 
world averages p^ . 
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TABLE III. Current world average and measurements for T^i from r — > KttttVt. 



^^-1(1270); GeV 1^^-1(1400); GeV 



PDG 1 


14 


1 


0.09 ±0.02 


0.174 ±0.013 


ALEPH H] 
DELPHI 11 
this analysis 


0.37 ±0.10 
0.19 ±0.07 
0.26t|]:[j? ± 0.08 


0.63 ±0.12 
0.31 ±0.08 
0.30l|];?j'±0.14 



B. Values of /1270 and fp 

As calculated in section VII, our data indicate that there is slightly more Xi(1270) than 
i^i(1400) in the axial vector current of the r~ K^tt^tt^Ut decay (/1270 = 0.66 ± 0.19 ± 0.13). 
Other experiments have also investigated the relative contributions of the two Ki resonances to 
r — suVr- One of the first measurements of the r — > KiVr branching fractions was performed by 
the TPC/2-i collaboration in 1994 ||]. Their results are B{t i^i(1270)z^) = 0.4lt[]:^^% and 
B{t i^i(1400)z^) = 0.76t{];3^%, giving the fraction /1270 = 0.35l[j:^^. The results of the TPC/2-f 
experiment suggest that the decay proceeds mostly through i^i(1400) although their errors are too 



large to draw firm conclusions. The latest branching fraction measurements by CLEO [17| and 
ALEPH [|^J2| as well as this analysis suggest i^i(1270) dominance. An analysis of the K~'k^ and 
7r~7r+ substructure in r~ K~Tr+TT~i^r allowed ALEPH [|l9| to determine /1270 = 0.41±0.19±0.15, 
based on the known branching fractions of the Ki resonances to K*tt and Kp. Recent measurements 
therefore favor i^i(1270) dominance in t~ K^tt^tt^Ut-. 

Calculating the amount of p in Ktttt from the fit parameters we find fp = 0.48±0.14, close to the 
measurement by ALEPH [l^ of fp = 0.35 ± 0.11. This number also agrees with the measurement 
by ALEPH of another related decay channel, r~ — > K^tt' where the component K p in the 
intermediate state is found to be (64 ± 9 ± 10)%, approximately twice that of K~ (P as expected by 
isospin symmetries Q. 



C. Ka — Kb mixing 

From our result for the ratio of r — > KiI't decay amplitudes, information about the mixing 
of the Ka and Kf, eigenstates can be derived. The mixing between Ka and Kj, is traditionally 
parameterized in the following way 

(1400) = Ka cos Ok - Kb sin Ok 
i^i (1270) = Ka sin Ok + Kb cos Ok 

In the case of exact SU{3)f symmetry, the second-class current r KbVr is forbidden and only 
Ka is produced. However, due to the difference between the masses of the up and strange quarks 
we may expect symmetry breaking effects of order \5\ = (m^ — m^jyplirns ± m^) ~ 0.18. Then, 
instead of pure Ka a linear combination \Ka) — b\Kb) is produced and the ratio of decay rates of 
the K\ resonances can be written as 

B{T^Kx(m^)v) 

b{t Ki(im)iy) 



sin Ok — cos Ok 



cos Ok + S sin Ok 



X 



(16) 



17 



In this expression, $ is the ratio of appropriate kinematical and phase space terms and is calculated 
by numerical integration. With the parameters measured in this analysis the ratio of branching 
fractions (Eq. |l^ is written as 

B{t ^ Ki{mQ)vr) ^ A^Ra + C^Rc , . 

From Eqs. ([T^-(p^), solutions for Ok can easily be found: 

(a) Ok = (69 ± 16 ± 19)° for 5 = 0.18, 

(b) Ok = (49 ± 16 ± 19)° for 5 = -0.18. 

There is a second pair of solutions that has opposite sign and the same magnitude. 

One can also calculate 6k using the current experimental information on the masses and branch- 
ing fractions of ii'i(1270) and i^i(1400), independent of their production in r-decays. There are 
two possible solutions, 6k ~ 33° and 6k ~ 57° Our result has the same two- fold ambiguity 
and is consistent with this calculation. 



X. SUMMARY AND CONCLUSIONS 

In this analysis we have measured the relative fractions and parameters of the Ki resonances 
in T~ — > K~TT'^'K~Vr dccays. These measurements are made within the framework of the model 
described in Sec. 0. Briefly, we assume that ifi(1270) and iCi(1400) saturate the if~7r+7r~ 
spectrum, and consider only the interference inherent in the Breit-Wigner mass distributions in 
calculating the relative t~ i^i(1270)z^T- and — > i^i (1400)1^,- branching fractions. Our param- 
eterization of the axial vector form-factors is different from Q in two respects - our form-factors 
are motivated by isospin relations, and we assume the Wess-Zumino anomaly to be negligible, as 
described in Sec. 0. We find /1270 = 0.66 ± 0.19 ± 0.13 and fp = 0.48 ± 0.14 ± 0.10, with /1270 
and fp defined as the -ftri(1270) and Kp fractions in r — > K~'K^Ti^Ur. 



These measurements agree well with the recent results from CLEO and ALEPH (see Sec VII). 
Our data slightly favor i^i(1270) dominance in production of the i^~7r^7r^ final state. The widths 
that we extract for the Ki resonances are considerably larger than previously tabulated values 



|14|. We also calculate the Ka — Kh mixing angle, finding 6k to be consistent with theoretical 



expectations. 
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